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Abstract 

We apply the method of QCD sum rules in the presence of the external electro- 
magnetic F^ v field to the problem of the electromagnetic decay of the various 
vector mesons, such as p ^ 7r + 7, K* — » K + 7, and rf — * p + 7. The induced 
condensates obtained previously from the study of baryon magnetic moments 
are adopted, thereby ensuring the parameter-free nature of the present calcula- 
tion. Further consistency is reinforced by invoking the various QCD sum rules 
for the meson masses. The numerical results on the various radiative decays 
agree very well with the experimental data. 
PACS Indices: 12.38.Lg; 13.40.Hq; 14.40.Aq 



Quantum chromodynamics (QCD) is believed to be the underlying theory of strong 
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interactions. However, little is known concerning the structure of a meson (for which 
the naive quark- ant iquark picture is at most only the leading approximation), since QCD 
is nonperturbative at the typical hadronic scale. In the absence of the first principle 
calculation, the electromagnetic decay of a vector meson, such as p — > 717, may be treated 
via an effective Lagrangian [[l], ||. In the vector meson dominance model (VDM) [0, 
an effective Lagrangian with small flavor 577(3) symmetry breaking was used while in 
chiral perturbation theory the Bardeen-subtracted Wess-Zumino action, together with 
symmetry-breaking corrections Fx/F n m F^jF^ « 1.2 and a given u-<j) mixing parameter 
e, was employed to predict radiative vector-meson decays. As pointed out in 0, VDM 
may work to the level of within 20% for the flavor SU(2) sector but is expected to be 
considerably worse for flavor SU(3). Another approach, which is the focus of the present 
paper, is made possible via the method of QCD sum rules || in which the nonperturbative 
QCD physics is incorporated systematically as power corrections in the short- distance 
operator product expansion (OPE). Results from QCD sum rules will certainly serve as 
another useful reference point concerning radiative decays of vector mesons. 

To be specific, we consider the electromagnetic decay of a vector meson in the flavor 
577(3) sector, i.e. in the sector of (u, d, s). Using p — > it + 7 as the illustrative example, 
we may describe the decay process by the effective Lagrangian, 



with Pcm is the pion momentum in the CM frame (i.e. in the rest frame of the vector 
meson) . 

In the original work of Shifman, Vainshtein, and Zakharov 0, the nonperturbative 
QCD effects are incorporated through the various non-zero condensates in the nontrivial 
QCD vacuum. Ioffe and Smilga [|J and, independently, Balitsky and Yung || studied 



C^^e^e^A^p^n, 



(1) 



where is the electromagnetic field. The corresponding partial width is 




(2) 
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baryon magnetic moments by extending the method of QCD sum rules in the presence 
of the external electromagnetic field. It turns out that the QCD sum rule method as 
obtained in jj| or can readily be applied to the problem of the electromagnetic decay 
of the various vector mesons, such as p — > n + 7, K* — > K + 7, and 7/ — ► p + 7. As 
a result, we wish to pursue along this direction a little further, hoping to obtain results 
which may help to clarify certain issues in relation to the effective Lagrangian approach 

We recall J§ g] that, in the method of QCD sum rules, the two-point correlation 
function in the presence of the constant electromagnetic field may be introduced 

as follows, 

IL» = 1 J rf 4 xe^(0|T[j PM (x), A (0)]|0) F , (3) 
where we introduce the interpolating fields, 

j Pfl (x) = u a (x)-i lJL d a (x) , j n (x) = u a (x)i>y 5 d a (x) , (4) 

with a the color index. As in the QCD sum rule analyses for nucleon properties, we may 
introduce the overlap amplitude (coupling) of the meson current to the (physical) meson: 

(0|jp» = , (5) 

(0|j.|vr) = A 7r , (6) 

where is the polarization 4— vector of the p meson. 

Embedding the system in the constant F^ u field and introducing intermediate states 
we express the correlation function at the hadronic level as follows, 

n» = aa^i - 1 - - 1 2 6^ ct f°V + nf-^(p) + • • • . (7) 

2m p p z — nip p z — p 
The prey coupling is defined through the relation: 



(K l ( P 'M m K(p)) = ^ lm e^e x ^ Puqa K( q 2 ), (8) 
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where K(q 2 ) is the form factor normalized such that K(0) = 1. The superscripts I and m 
are isospin indices. Note that, in Eq. (7), we write down explicitly the leading p — > 717 
contribution, label the next transition term by Ti p *~ i ' K1 {p) 1 and denote the continuum 
contribution simply by the ellipsis. 

The external field may induce changes in the physical vacuum and modify the 
propagation of quarks. Up to dimension six (d < 6), we may introduce three induced 
condensates as follows. 

(0|ga^g|0)^ = e g e X F^(0\qq\0), (9) 

9c(0\q Y GU\°)^ = e g eKF^(0\qq\0), (10) 
yn 

9c^x«{mi^Gl a q\U) Fllu = ie q eCF^(0\qq\0) , (11) 

where q refers mainly to u or d (with suitable modifications for the s quark), e is the 
unit charge, e u — |, and = — |. In studying baryon magnetic moments, Ioffe and 
Smilga obtained the nucleon anomalous magnetic moments fi p = 3.0 and fi n = —2.0 
(±10%) with the quark condensate susceptibilities x ~ — 8 GeV 2 and k, £ quite small. 
Balitsky and Yung || adopted the one-pole approximation to estimate the susceptibilities 
and obtained, 

x = -3.3 GeV" 2 , k = 0.22, £ = -0.44 . (12) 

Subsequently Belyaev and Kogan |J extended the calculation and obtained an improved 
estimate x — — 5.7GeV~ 2 using the two-pole approximation. Chiu et al [0] also estimated 
the susceptibilities with the two-pole model and obtained 

x = _4.4GeV~ 2 , k = 0A, £ = -0.8. (13) 

Furthermore, Chiu and co-workers |7j re-analysed the various sum rules by treating Xi 
k, and £ as free parameters to provide an overall fit to the observed baryon magnetic 
moments. The opitimal values which they obtained are given by 

x = -3GeV~ 2 , k = 0.75, £ = -1.5. (14) 
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We note that in these analyses the susceptibility values are consistent with one another 
except the earliest result x = ~ 8 GeV -2 in which is considerably larger (in magnitude). 
In what follows, we shall adopt the condensate parameters x = — 3.5 GeV -2 , 2k + £ m 
with k, £ quite small, which represent the average of the latter values discussed above. 

On the other hand, the correlation function ^(p) may also be evaluated at the quark 
level. 

IL» = -if d A x^Tr[iSl b {x) lfl iS b d a {-x)i l5 l (15) 

where iS ab (x) is the coordinate-space propagator in the presence of the field [|], 

Equating the correlation function obtained both at the quark level (l.h.s.) and at the 
hadron level (r.h.s.), performing the Borel transform on both the sides (to accentuate the 
contribution from the leading p — > tt-/ contribution), we finally obtain the QCD sum rule 

for g piT1 : 

(e u + e d )a q ^ [-3 X + (2 + 2k + 0^] + • • • 

ml B $ (16) 
= f^4^(e ^-e ^) + n p ,_ 7 (M), 

where a q = — (27r) 2 (0|gg|0) and Mb is the Borel mass (with M B playing essentially the 
role of p 2 ). Note that in this sum rule (16) the dominant contribution comes from the 
induced quark condensate characterized by the susceptibility x- The contribution from 
p* — > 7T7 is suppressed by the factor (^-) 3 and so is negligible to within 10% (which is 
the typical accuracy of the QCD sum rule approach). The sum rule ( |I6"D is stable within 
the working interval M B ~ m 2 . 

We proceed to note that in the framework of QCD sum rules the overlap amplitude 
A p can be determined in a self-consistent manner by making use of the mass sum rules 
for the p and tt mesons 0. |[. 



where E\ = 1 — (1+ jrr)e (with S% = 1.5GeV 2 ) is the factor arising from a quark-level 

B 

approximation for the continuum, and b = (27r) 2 (0|^-G^ (g G™ (3 |0) is the gluon condensate. 
It was pointed out in |H| that the instanton contributions may invalidate the usual sum 



rule techniques for the pseudoscalar current. Nevertheless, it was suggested in [||, |8| 
that Eq. (18) may still provide a rough estimate of A„- with the experimental pion mass 
as input (despite the fact that it cannot be used to predict the pion mass). Moreover 
there are many careful analyses of this correlator which yield very reasonable results for 
A, 0, U, H. We therefore adopt the value from the QCD sum rule for A„- with the 



physical mass as the input. We obtain the estimates A^ ~ 0.17 ± 0.03 GeV 2 and 



A p ~ 0.17 GeV 2 . An independent phenomenological analysis in [[14]] yields A p = / ( 



0.17GeV 2 , A,r = A — m + m = 0.20±0.04GeV 2 , values consistent with those from the present 
QCD sum rule approach. 

Using our result on X p and A„-, the observed masses, and the condensate parameters 
discussed earlier, we finally obtain from the QCD sum rule (16), 

^ 7 = 0.59, (19) 

which is very close to the value g pnl = 0.58 as may be extracted from the partial width 
for the process p — > ivy. 

It is clear that the method described so far for p — > 777 can be generalized in a 
straightforward manner to other electromagnetic decays of vector mesons, among which 
we have considered additional processes as listed in Table I. The working interval for the 
Borel mass for the corresponding sum rules is M\ ~ m 2 ar , where m par is the mass of the 
parent particle in the decay process. In the case for the meson containing a strange quark 
or antiquark, we have taken into account the contributions due to the nonzero quark mass 
m s 150MeV, which in the cases such as K* — > Kj may give rise to corrections of order 
of about 20%. We also use (0|ss|0) ~ 0.8 x (0|mm|0) (as suggested by chiral perturbation 
theory) which is essential for the K* — > Ky channel. 
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In Table I, we present the experimental values for the various electromagnetic decays 
of vector mesons, the predictions from VDM [I], those from the effective Lagrangian 
approach 0, and what we have obtained from the QCD sum rules (this work). The 
results shown in this table indicate that the QCD sum rules yield predictions which are 
in better overall agreement with the experimantal values. 



TABLE I. Comparison between experimental data H] and predictions from the Vector 
Meson Dominance Model (VDM) effective chiral Lagrangian approach M, and QCD 
sum rules. The unit is KeV. 



u> 

P 
P 



K°* 
K+* 

V' 

i 



7T7 

• 777 
7T7 
777 

VI 

ii 

> fsT°7 

> if+7 

* PI 
■ CJ7 

7T7 



experimental 
value 



716 
7 

68 
58 
57 
< 1.S 
116 
50 
60 
6 

5.8 



Vector Meson 
Dominance 



888.0 
5.4 
63.0 
55.0 
71.0 
0.24 
149.0 
68.0 
119.0 
12.5 
6.1 



effective chiral 
Lagrangian 



800 
5 

80 
38 
68 
1 

117 
29 
77 
7 
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QCD sum 
rules 



720 
7.1 
65 
48 
80 
0.67 
102 
60 
60 a 
6.1 
5.8 b 



a Used as the input to determine rj-7]' mixing angle 9. 

b Used as the input to determine cf)-uj mixing parameter e. 
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We note that all parameters have been determined consistently in the framework 
of QCD sum rules as illustrated before. In particular, the various overlap amplitudes 
(couplings) are determined from the various mass sum rules: A^ = 0.17 ± 0.03GeV~ 2 , 
\ K = 0.23 ± 0.03GeV~ 2 , \ p = 0.17 ± 0.01GeV~ 2 , \ K * = 0.21 ± 0.02GeV~ 2 . The above 
predictions and the other overlap amplitudes given in the following after taking into 



account the mixing effects are in good agreement with those in [pL4 ] . Note that the value 
for A^ is of some specific interest in view of the £7a(1) problem. 
In the case of the <p-uo mixing, we introduce 



and 



where 



\uj) = cose|co> ) — sine|0 o )) 10) = shield) + cose|0 o ) , (20) 



j u — cos e - sin e , j<f> = sin e + cos e . (21) 



\co ) = ^(\uu) + \dd)), \<j> ) = \ss). (22) 

Using the interpolating field in (21), we derive the QCD sum rules for the relevant radiative 
decays and the mass sum rules for u and 4> mesons. We then use the observed value for 
the width of the decay <fi — > tc-/ to determine the mixing angle. We obtain the mixing 
parameter e = 0.063 ±0.005 and the overlap amplitudes A w = 0.16±0.01GeV~ 2 and A^ = 
0.24±0.02GeV" 2 . As the Borel mass is varied in the interval 0.6GeV 2 < M 2 B < l.lGeV 2 , 
the mixing parameter e changes from 0.068 to 0.058. Note that, as a consistency check, 
our result is consistent with the value e = 0.079 as determined from cf) — » pn 0. 

Furthermore, we use the experimental partial width for the decay process r( — > p7 
as the input in the derived QCD sum rules to extract the r] — ^ 8 mixing angle, which is 



parametrized as in JT5 



\rj) = cos 9\i] 8 ) - sin 9\rj ), \rj') = sm0\r]g) + cosd\rj ) , (23) 
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and correspondingly for the interpolating fields, 

j n = cos 6j m - sin 6j vo , j n > = sin 8j m + cos 6j vo . (24) 

At the same time we obtain the 77, rj' mass sum rules with the interpolating fields in (24) 
and use the physical meson masses to calculate the overlap amplitudes. We then use the 
extracted mixing angle to predict other processes. Our results are = 0.23±0.03GeV~ 2 , 
A^/ = 0.33 ± 0.05GeV -2 . 9 = —19° ± 2°, a value consistent with that extracted from the 



decays rj — ► 77 and rf — > 77 [ 15 1 and in || . We find that such value on the 770 — r)$ mixing 
is essential for the overall agreement with the experimental data. 

We wish to emphasize the parameter-free nature of the present calculation, in that all 
parameters are determined consistently in the framework of QCD sum rules, with a single 
set of the values for the induced condensates obtained previously from studies of baryon 
magnetic moments. We note that the uncertainty of the overlap amplitudes (couplings) 
derived in this way might be as large as 10% in some cases and the method of QCD 
sum rules, being an operator product expansion used for moderate Q 2 , may have some 
inherent error of order of 10%. With this in mind, we may conclude that our results as 
shown in Table I agree very well with the experimental data. 

This work was supported in part by the National Natural Science Foundation of China 
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